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ABSTRACT 


We calculate the abundances of electrons and ions in the hot ( ^ 500 K), dusty 
parts of protoplanetary disks, treating for the first time the effects of thermionic 
and ion emission from the dust grains. High-temperature ionization modeling has 
involved simply assuming that alkali elements such as potassium occur as gas- 
phase atoms and are collisionally ionized following the Saha equation. We show 
that the Saha equation often does not hold, because free charges are produced 
by thermionic and ion emission and destroyed when they stick to grain surfaces. 
This means the ionization state depends not on the hrst ionization potential of 
the alkali atoms, but rather on the grains’ work functions. The charged species’ 
abundances typically rise abruptly above about 800 K, with little qualitative de¬ 
pendence on the work function, gas density, or dust-to-gas mass ratio. Applying 
our results, we hnd that protoplanetary disks’ dead zone, where high diffusivities 
stifle magnetorotational turbulence, has its inner edge located where the tem¬ 
perature exceeds a threshold value ~ 1000 K. The threshold is set by ambipolar 
diffusion except at the highest densities, where it is set by Ohmic resistivity. We 
hnd that the disk gas can be diffusively loaded onto the stellar magnetosphere 
at temperatures below a similar threshold. We investigate whether the “short- 
circuit” instability of current sheets can operate in disks and hnd that it cannot, 
or works only in a narrow range of conditions; it appears not to be the chondrule 
formation mechanism. We also suggest that thermionic emission is important for 
determining the rate of Ohmic heating in hot Jupiters. 

Subject headings: accretion, accretion disks, planetary systems: protoplane¬ 
tary disks, solar system: formation 
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1. Introduction 

Protoplanetary disks’ evolution is governed by the extraction of orbital angular 
momentum needed for accretion to proceed. Magnetic forces are among the prime 
candidates to drive the extraction. Modeling the disks therefore requires knowing how 
readily the magnetic fields diffuse through the material. Hence the ionization state of the 
gas is key. 

In particular, the magnetorotational instability (MRI) is suppressed unless the electron 
fraction is Xe = ■ne/'nH 2 ~ 10“^^ (Jin 1996; Gammie 1996), where nH 2 is the number density 
of H 2 molecules. Launching disk winds requires ^ 10“^^ — 10“^° (Bai & Stone 2013). 
The gas ionization likewise bears on the existence of lightning (Desch & Cuzzi 2000) and 
current sheets (McNally et al. 2013), two proposed mechanisms for transiently heating 
and melting chondrules (millimeter-sized inclusions in chondritic meteorites) in the solar 
nebula. Assessing the extent of the MRI and the heating of chondrules requires knowing 
the ionization rates, especially in the disks’ denser midplane regions. 

Ionization in the cooler parts of disks, where temperatures T < 500 K, has been 
studied in detail. Several non-thermal processes are important at these low temperatures, 
with ionization rates per volume often expressed as C^H 2 - Far-ultraviolet (FUV) photons 
from the central star can ionize gas at a rate Cuv ~ 10“® s“^ at 1 AU, but are attenuated 
by column densities ^ 0.1 gcm“^ (Perez-Becker & Chiang 2011). Similarly, X-rays from 
the star can ionize at rates Cx ~ 10“^^ s“^ at 1 AU, but are attenuated by column densities 
< 10 gcm“^ (Glassgold et al. 1997; Igea & Glassgold 1999; Ercolano & Glassgold 2013). 
Galactic cosmic rays, if not shielded by stellar winds (Gleeves et al. 2013a), can ionize at 
rates Cgcr 10 ^'^s but are attenuated by column densities ~ 100 gem ^ (Umebayashi 
& Nakano 1981). In the innermost few AU of a protoplanetary disk, the expected column 
densities are 3> lO^gcm”^, so that the dominant non-thermal ionization sources are /3 
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particles from radioactive decays of species like leading to Crad ~ 10“^^ s“^ (Umebayashi 
& Nakano 1981); or ^®A1 (if present), leading to Crad ~ 10“^® s“^ (Cleeves et al. 2013b). In 
the disks’ cold midplane regions, this is the expected ionization rate, and it is quite low. 

In chemical equilibrium in the dense midplane gas, the electron density is found by 
balancing the ionization rate against adsorption, and subsequent recombination, of electrons 
on grain surfaces. Assuming monodisperse grains with radius Ogr = 1 /rm, internal density 
Ps = 3gcm“^, and dust-to-gas mass ratio Pgr/Pg = 0.01, the equation describing this 
balance is 

2 /SkTV''^ 

CnHa = TlgrTra^r^e I - I *Se, (1) 

yTTTTT-e / 

(Sano et al. 2000), where ngr = (pgr/Pg) (l-4nH2^^H2) /i^'^Ps^gr/^) is the number density of 
grains, S'e the sticking coefficient of electrons colliding with grains, and other variables have 
their usual meanings. From this we hnd Ug 10 ^ cm ^ (assuming^ = 10 ^®s T = 300 K 
and S'e = 0.6). For a typical midplane gas density inside a few AU, ^ 10^"^cm“^ 
(Weidenschilling 1977), Xg ^ 10“^^, a value too low to be of signihcance to the MRI. 
Recognition of this fact led Gammie (1996) to suggest that protoplanetary disks’ interiors, 
especially the parts within a few AU of their stars, are “dead zones” where the MRI cannot 
act. 


Jin (1996) and Gammie (1996) were the first to recognize that very close to the central 
star, protoplanetary disks may be hot enough to sustain thermal ionizations, recoupling 
the magnetic held to the gas. Specihcally, they pointed out that collisional ionization of 
gas-phase potassium atoms at temperatures > 1000 K would generate free electrons and 
ions. The threshold temperature of 1000 K is based on calculating the abundances of K"*" 
ions and electrons using the Saha equation: 

nK+ g+ 2 f 2'KmgkT \ / IP \ 

^ ^ ^ ^ V ■ 

Here = 1 and go = 2 are the statistical weights of the (ground states) of K+ and K°, k 
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the Boltzmann constant, h Planck’s constant, rrie the electron mass, and IP = 4.34 eV the 
hrst ionization potential energy of potassium. Assuming nK+ = Ue and that the abundance 
of K relative to H 2 in a solar-composition gas is xk = 3.0 x 10“^ (Lodders 2003), half 
of the K atoms are ionized when T = 1715 K, for a density UHa = lO^^cm”^ likely to 
occur near the midplane. More significantly for the MRI, Xe = 10 when T = 865 K, 
meaning substantial ionization would be expected if potassium were in the gas phase at this 
temperature. According to Lodders (2003), the 50% condensation temperature of potassium 
is 1006 K (at P = lO^'^bar, or = 6 x lO^^cm”^; the pressure dependence is weak). So it 
is fairer to say that above 1000 K the Saha equation predicts signihcant thermal ionization 
of potassium, and for temperatures between about 850 and 1000 K, thermal ionization may 
or may not lead to signihcant ionization. 

Other species contribute less strongly to thermal ionization than potassium. Because 
of the exponential dependence on the hrst ionization potential (IP), one might expect other 
alkali atoms to contribute to thermal ionization. The hrst ionization potentials of Li, Na, K, 
Rb, and Cs are, respectively, 5.39, 5.14, 4.34, 4.18, and 3.89 eV. These elements’ abundances 
relative to H 2 are, respectively, 4.57 x 10“®, 4.73 x 10“®, 3.04 x 10“'^, 5.51 x 10“^°, and 
3.02 X 10 (Lodders 2003). The hrst ionization potentials (PIP) and cosmochemical 
abundances of many elements are depicted in Figure 1. Repeating the analysis above, the 
temperature at which each alkali element contributes enough electrons to yield Xe = 10“^^ 
(at nu 2 = lO^^cm”^) is 1149, 975, 865, 932, and 852 K, respectively. Considering that 
cesium is in the gas phase above 800 K whereas potassium is not (Lodders 2003), one might 
expect Cs to be the dominant contributor of ions and electrons at temperatures between 
850 and 1000 K; however, the Cs is 10^ times less abundant, so that even if all Cs were 
ionized and only 0.1% of all K atoms were in the gas phase, potassium would still dominate 
over cesium. Thermal ionization of potassium is likely to be more important than that of 
other alkalis, and potassium is in any case an excellent proxy for other alkalis. Throughout 



this paper we will consider the case of potassium to represent the contributions of all the 
alkali elements to thermal ionization. 



7 


> 

(U 


cc 

■p 

c 

(U 

O 

Dh 

c 

o 

■l-l 

cC 

N 

• fH 

c 

o 


m 


15 


“I-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-r 


10 


® _ 

€) » ei 


(3 


XO ® 


A 




Cs 




X* 

Rb 


) ()< o 


x-O 


Meteoritic_Photospheric 




CSC 

N x-^ ® g, «l 

^ ca o< 


® ® ® 

^Al 

•< 

» Na 
K 


'Mg 


0 2 4 6 8 10 12 

Solar Abundance by log(Number), H=12 


Fig. 1.— Cosmochemical abundances of various elements and their first ionization potentials 
(FIPs). Abundances (on a scale where the number of H atoms is 10^^) are derived by Asplund 
et ah (2005) from analyses of carbonaceous chondrites (“meteoritic”) or from spectrometry 
of the Sun (“photospheric”). Producing free electrons by thermal ionization requires the 
presence of elements with high abundance and low FIP, i.e., those elements toward the 
plot’s lower right corner. The exponential sensitivity of the ionization fraction to FIP means 
potassium (K) is most favored to provide free electrons by thermal ionization, with Na, Rb, 
and Cs potentially also important—if these elements occur in the gas phase rather than 
locked up in minerals. 
















A criticism of this thermal ionization picture is that the Saha equation is valid only in 
local thermodynamic equilibrium, or LTE (Desch 1998). It relies on the collisional ionization 
of the gas-phase K atoms, e.g., H 2 -|- K° —)■ H 2 -|- K+ -|- e“, being in detailed balance with 
its inverse reaction, three-body gas-phase recombination, H 2 -|- K"*" -|- e“ —)■ H 2 -|- K°. But 
under conditions typical for protoplanetary disks, electrons are more likely to recombine 
with ions on the surfaces of grains. Desch (1998) argued that the electron density must be 
set by balancing collisional ionizations against collisions of electrons with dust grains. Once 
adsorbed on the grain surfaces, electrons and ions quantum tunnel from one adsorption 
site to the next, wandering across the surface until they hud each other and recombine. 
The resulting neutral K atoms can then escape the grain. Dust grain surface-catalyzed 
recombinations are much faster than gas-phase recombinations, as the frequency of 
colliding with dust grains is ~ ngrTTOg,. (8/cT/7rme)^'^^ S'e ~ 0.1s“^ (at T = 1000 K and 
nHa = 10^"^cm“^, and assuming grains with radius 1 /rm and a dust-to-gas mass ratio 0.01), 
whereas the frequency of recombining with an ion is ~ ^k+ ~ 10“^^ s“^, where A ;_2 

is the recombination coefficient (see below), and we have assumed the ion density matches 
the electron density, and = 10“^^. Not until ^ 10“^ do gas-phase recombinations 
dominate over grain-catalyzed recombinations. The dominance of grain effects invalidates 
the use of the Saha equation for potassium as a means of calculating the electron density. 

The dominance of grain processes also points to an additional ionization source, not 
previously considered in the context of protoplanetary disks. The process that is the inverse 
of electrons and ions colliding with grains, resulting in neutral atoms leaving the grain, is 
the emission of electrons or ions following collisions of neutral atoms with grains. These 
processes are more commonly called thermionic emission and ion emission. When solids are 
heated to temperatures > 500 K, they emit electrons and/or ions. The emission of electrons 
depends on the work function of the solids comprising the grains, which also affects whether 
evaporating atoms such as potassium are emitted as neutral atoms or ions. Including the 
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thermionic and ion emission completely changes the role of dust grains in the ionization of 
the gas: instead of mopping up free charges, they generate them. 

In this paper we examine the effects of thermionic and ion emission on the ionization 
state of high-temperature, dusty gas. In §2 we discuss the microphysical mechanisms of 
thermionic and ion emission, as well as the evaporation and condensation of K atoms, and 
provide quantitative rates and likely input parameters. In §3 we create a simple model 
treating evaporation and condensation of K atoms and ions, including ion emission and 
thermionic emission. We show that, because of the high work function of silicate grains, 
ion emission creates signhcant ionization at temperatures around 800 K, possibly as low as 
600 K. In §4 we discuss the impact on several aspects of protoplanetary disks. In particular, 
we address where the dead zone’s inner edge should he, what temperatures are needed for 
disk gas to diffuse into the star’s magnetosphere, and whether including this new ionization 
source can enable the “short-circuit instability” of McNally et ah (2013). Thermionic and 
ion emission profoundly affect the ionization state of high-temperature gas in disks. 


2. High-Temperature Chemical Network 
2.1. High-Temperature Ionization Effects 

In the dusty gas of protoplanetary disks, we cannot always use the Saha equation, and 
instead must check the rates of various ionization and recombination processes. What’s 
more, protoplanetary disks are marked by transient temperature changes, and we are 
interested in how the ionization responds to the changing conditions. Thus motivated, in 
this section we lay out several processes’ kinetic rates and construct a chemical network 
to hnd the steady-state abundances of electrons and ions as well as the timescales for the 
ionization chemistry to reach equilibrium. 



Our chemical network includes five species: electrons, molecular or atomic ions, K 
atoms bound to dust grains, gas-phase neutral K atoms, and gas-phase K ions. Energetic 
particles (cosmic rays, radioactive decay products, etc.), and at high temperatures also 
gas-phase reactions, ionize the gas molecules; ions and electrons recombine in the gas phase. 
We also include dust grains, calculating their steady-state charge. All the other species 
can collide with and condense on or chemically react with the grains. The grains can lose 
K atoms at high temperatnres, and can eject electrons and ions by the thermionic effect 
and by ion emission, respectively. Before discussing the chemical network, we review these 
chemical processes. 


2.2. Gas-Phase Collisional Ionization and Recombination 

As discussed above, the H 2 molecules’ ionization can come from energetic particles 
snch as Galactic cosmic rays. X-ray photons, or particles from radioactive decays. The 
nnmber of ionizations per volnme per time is defined as ( UHa • We are most interested in 
protoplanetary disks’ interiors, so we assume that the cosmic rays and X-rays are shielded, 
leaving only ionization by the decay of radioactive elements. If short-lived radionuclides 
are not present, then decay of as well as dominate, yielding 

C ^ 6.9 X 10“^^ s“^ (Umebayashi & Nakano 1981). If the short-lived {ti /2 = 0.71 Myr) 
radionnclide ^®A1 is present at levels inferred from meteorites (MacPherson et ah 1995), 
then C, ~ 10“^®s“^ (Umebayashi & Nakano 2009). Onr focus here is on high-temperatnre 
ionization, so we arbitrarily take C, = 6.9 x 10“^^ s“^ in what follows. (If ionizations by ^®A1 
are significant, the ion and electron densities wonld rise sharply above 800 K instead of 
700 K.) The ions qnickly exchange charge with other molecules to produce ions snch 
as HCO"*", which then also exchange charge with free atoms, creating ions snch as Mg^. 
We assnme these steps are faster than the other reactions, so each ionization of H 2 qnickly 
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yields an atomic ion. We do not consider K"*" as one of these ions, because other species 
like Mg are more abundant. Instead we assume K"*" ions are produced mainly by thermal 
ionization. 

Atomic ions and free electrons can recombine in the gas phase. We assume a rate 
per volume where f3 = 3 x 10“^^ cm^ s“^, T measured in K (Oppenheimer & 

Dalgarno 1974). Ions and free electrons can also be adsorbed onto grain surfaces, at rates 

77i,ads = ningrTTOg, {SkT/nniiY^'^ Si Ji = Ui z/i,ads, (3) 

and 

77e,ads = {3kT/ Si Je = ne Z/e,ads, (4) 

where J, and Je reflect modihcations to the collision cross section due to the grain charge 
(Draine & Sutin 1987). 

Thermal ionization of alkali atoms occurs if they are in the gas phase and suffer 
collisions with ambient molecules, mostly H 2 . This process can be represented by the 
reaction 

H2 + K°^H2 + K+ + e-. 

In the gas phase, electrons and ions can radiatively recombine, as follows, 

K+ + e" ^ K° + 7 , 

but in the regions of protoplanetary disks where temperatures are > 500 K (the midplane 
regions well inside 1 AU), the gas densities are likely high enough that three-body 
recombinations dominate: 

H 2 + K+ + e- -> H 2 + K°. 


We justify this statement below. 
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Literature on the collisional ionization and recombination rates is sparse. Pneumann & 
Mitchell (1965) present a kinetic rate, but cite no measurements or sources. Their rate also 
appears to conflict with that provided by Ashton & Hayhurst (1973). Since the latter is 
experimentally determined by analyzing the chemical reactions in flames, it is our preferred 
value. Following Ashton & Hayhurst (1973), collisional ionization takes place at a rate per 
unit volume 


77gas,ion ^2 ^H 2 j 


where 


k2 


(9.9 ±2.7) X exp 



—1 

cm s 


(5) 

( 6 ) 


and T is measured in K. Here IP = 4.34 eV for K, and IP/A; = 50350 K. The 3-body 
gas-phase recombination rate is 


77gas,rec k—2 "^^2 ^K+ 


(7) 


where 

k -2 = (4.4 ± 1.1) X 10"^^ cm^^ s"\ 
and again T is measured in Kelvin (Ashton & Hayhurst 1973). 


( 8 ) 


The two coefficients k 2 and k _2 must be related to each other. In the absence of dust, 
these would be the only rates affecting the ionization fraction, and should be in detailed 
balance with each other in local thermodynamic equilibrium. The ionization fractions found 
using these coefficients thus should conform to the Saha equation. Hence, 77gas.rpn = 
and nK+ Ue /i ;_2 = k 2 , and therefore we expect 


k_ 


3/2 


k2 UK+Ue 2 (/+ /'27rmekT\ 

j exp 




-2 ^Ko go 
The values from Ashton & Hayhurst (1973) yield 


= 2.415 X IQi'T^A exp (9) 


± = (2,3±0,8)xl0-rV^exp(-P). 


( 10 ) 
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consistent with the Saha equation. 

Recombinations by 3-body reactions proceed at a rate (/c _2 nna) ''^e, which is 
to be compared to the rate of radiative 2-body recombinations, /Inx+ne. For typical 
radiative recombinations, /d ~ 3 x 10“^^ cm^ s“^ (Oppenheimer & Dalgarno 1974). At 
T = 500 K, /c _2 ~ 10“^®cm®s“^, and if nHj = lO^^cm”^, then A;_ 2 nH 2 ~ 10“^^cm^s“^. For 
densities > lO^^cm”^, radiative recombinations of K"*" and electrons are less signihcant 
than 3-body recombinations, but we include both rates in our calculations. 


2.3. Thermionic Emission 


Thermionic emission, the ejection of electrons from heated solids, is a long recognized 
effect and the basis for the cathode ray tube. The rate at which electrons are emitted, per 
surface area of solid, is described by Richardson’s law, 

/ W\ 

( 11 ) 


, 47rme(A:T)2 f W 
]{T) = Ar --exp ( - 




kT 


where Ar is a dimensionless number of order unity depending on the solid material, and is 
typically ^1/2 (Crowell et ah 1965). 


Thermionic emission is sensitive to W, the work function of the material - that is, 
the energy needed to make an electron escape the solid. As the grains lose electrons and 
develop a positive charge +Ze, removing more electrons requires extra work to overcome 
the Coulomb force. The effective work function is 

IFeff = IF + —. (12) 

Ogr 

For conductors, the escape takes place from the top of the conduction band, and the work 
function is easily measured and commonly tabulated. For example, the work function of 
graphite is 4.62 eV (Jain & Krishnan 1952), Fe has a work function of 4.31 eV, and that of 
Ni is 4.50 eV (Fomenko 1966). 
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Unfortunately, the grains likely to occur in the solar nebula and other protoplanetary 
disks contain many insulating materials, making determining the work function difficult. 
The best analogs to the grains in the solar nebula, and presumably in other protoplanetary 
disks as well, are the unaltered (pre-accretionary) matrix grains in chondrites. These grains 
are composed of crystalline, Mg-rich olivines [Mg 2 Si 04 ] and pyroxenes [MgSiOg]; relatively 
Fe-rich amorphous silicates [(Mg, Fe) 2 Si 04 ] and [(Mg, Fe)Si 03 ]; and a smaller fraction (a 
few percent) of grains of metallic Fe,Ni and their sulhdes (Scott & Krot 2005). Metallic 
FeNi grains are conducting, with a work function presumably ~ 4.4 eV, but the silicates 
are insulators. Because of this, work functions of silicates are difficult to measure, and 
few data exist in the literature. Quartz (Si02) has a work function 5.0 eV, and mica (a 
phyllosilicate) has a work function 4.8 eV (Fomenko 1966), but the authors were unable to 
hnd any other direct measurements for specihc silicate minerals. Far better than quartz and 
mica as analogs for solar nebula dust, however, is the JSC-1 lunar regolith simulant (McKay 
et al. 1993), composed of a mixture of 40 wt% ferromagnesian olivine [(Fe, Mg) 2 Si 04 ], 

40 wt% plagioclase [NaAlSigOg and CaAl2Si208], the remainder comprised of Fe- and 
Ti-rich minerals like ilmenite (FeTiOs), anatase [Ti02], magnetite [Fe 304 ], hematite [Fe 203 ] 
and pseudobrookite [Fe 2 Ti 05 ]. Trigwell et al. (2009) examined how JSC-1 simulant is 
charged when brought into contact with different substances, and found it was negatively 
charged by Al (work function 4.28 eV), Cu (4.65 eV) and stainless steel (5.04 eV), but 
positively charged by the organic polymer PTFE (work function 5.75 eV). The acquired 
charge correlates with the work function of the other material, and the correlation strongly 
suggests JSC-1 itself has an effective work function somewhere between 5.0 and 5.4 eV. 
Feuerbacher et al. (1972) examined the photoelectric effect using JSC-1 simulant and found 
the threshold for photoelectron emission implies an effective work function close to 5.0 eV. 
Based on these measurements, we adopt a work function of 5.0 eV for solar nebula dust 


grams. 
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With the work function dehned, the rate of electron production per volume due to 
thermionic emission from dust is 

Weff 


.2 „ ._„2 X 47rme(/i;T)^ 


"^therm ^gr 47rClgj. J (7^) TZgp 47rClgj. 


exp 


(13) 


*3 V kT 

Again, for silicates we adopt W = 5.0 eV. Since data are lacking for the coefficient Ar, we 
adopt Ar = 1/2. 

The inverse of thermionic emission is electrons sticking on dust grains, which occurs at 
a rate Ue z^ads- In the absence of gas-phase processes affecting electrons, thermionic emission 
and its inverse could be expected to be in detailed balance, with T^therm = ads and 


. Anm.ek'^T'^ 

Ar -- exp 


'^e . I I 

4 \7rmeJ 


(14) 


V kT 

Note that since the electrons’ production and destruction rates both scale with the grain 
snrface area, the grain density and radii drop out of the balance equation. Solving for Ue 
yields 


Up = 


4Ar k 27rmekT\^^‘^ 

[~~h^ J 


2 JpS'e 


exp 


Weff 

' kT 


(15) 


This is identical in form to a Saha eqnation, involving a charged and nentral species X~^ 
and with statistical weights g~^ and g^, provided {nxo/nx+){g+/go) = ^r/{ 2Je Se) ~ 1. 
Potassium ion and electron densities essentially follow a Saha equation, but with an energy 
equal to Wgs, which in practice is typically in the range of 3 to 4 eV (see below), not the 
ionization potential IP = 4.34 eV. In this way the dust grains fnndamentally alter the gas 
ionization state, even in the high-temperature limit, where their effects are described by a 
Saha-like equation. 


2.4. Evaporation of Neutral and Ionized Potassium Atoms 

In addition to electrons, atoms can also emit ions. For the reasons outlined above, most 
relevant are K+ ions. At high temperatnres, potassinm exists mostly in the gas phase, with 
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K atoms that meet dust grains quickly vaporizing because of thermal vibrations. A single 
K atom bound to a grain surface vibrates on the lattice with a frequency v. We adopt a 
typical value u = 3.7 x appropriate for K bound to Pt (Hagstrom et al. 2000). In 

any vibration cycle, there is a probability of vaporization exp(—i7a//^T), where Eg, is the 
activation energy, equal to the binding energy, typically several eV. Dehning nK,cond as the 
number of K atoms in the grain minerals per volume of nebula, the vaporization rate per 
volume is 


TT-Kjevap ^KjCond ^ exp 



^K,cond ^evap‘ 


( 16 ) 


Balancing these vaporizations is the condensation of K atoms onto the grain surfaces. The 


total rate at which K atoms collide with dust grains is 77 ko,co11 + TT-K+.coih where 


/ \ 1/2 

"/^K+jColl '^K+ Ugi-TT Ogj, I I J\ = '^K+ ^K+,colh 

^ \TTmKj 


Ji being the focusing factor appropriate for ions (Draine & Sutin 1987) and tuk the mass of 
a potassium atom, and 


77ko,co11 Tlgr TT Ugj. 


f 8kT \ _ 

- = ^KO.coll- 

XTruiK/ 


( 18 ) 


The ratio of bound K atoms to gas-phase K atoms can be found be balancing these two 
rates: 


^KjCond 

riK 


~ Ugi. TTOgj. 


exp 

TTTUk / \ KI 


ipjpi) /skTy'’^ 


^'^dpsflgr/S yTrmKj 


o — i / I 

SiU exp (+ — 


(19) 


The fraction of potassium atoms in the gas phase (versus bound in grains) is strongly 
temperature dependent. 


At low temperatures, potassium is unlikely to be found in the gas phase because it 
will chemically react with silicates and be bound into minerals. In particular, potassium 
resides in feldspar (KAlSisOg) below its 50% condensation temperature ~ 1006 K (Lodders 
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2003). To mimic this behavior, we assume the binding energy of K atoms in silicates 
is Eg_ = 3.26 eV, so that nK,coii/''^K = 1/2 at T = 1006 K. Here we take a gas density 
UHa = 10^"^cm“^, a dust-to-gas ratio = 0.01, a grain density ps = 3gcm“^, a grain radius 
Ogr = 1 /im, a sticking coefficient S'! = 1, and neglecting the charge state of the K atoms. 


At temperatures above 1006 K, potassium atoms will evaporate from the surface, but 
whether they leave as neutral atoms or as ions depends on their ionization potential and the 
work function of the grains, according to the Saha-Langmuir equation. Among departing 
potassium atoms, the ratio of ions to neutrals obeys 


nK+ 9+ / fTefr-IP\ 

-= — exp +-—- . 

uro go \ kT J 

The fraction of leaving K atoms that are positive ions is then 


( 20 ) 


(wk+Zuko) 

1 + {ri^+/nKoy 


( 21 ) 


The rate at which ions evaporate from grains (per volume of nebula) is therefore 77K,evap/+, 
and the rate at which neutral K atoms do so is 77K,evap(l — /+)• 


2.5. Chemical Network 

We now construct a chemical network to calculate the densities of the hve species, 

Ue (free electrons), (gas-phase atomic ions), uk.coII (grain-adsorbed K atoms), nK+ 
(gas-phase K ions), and Uro (gas-phase neutral K atoms), as well as the charge Z on each 
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dust grain. The time derivatives of these quantities are as follows. 


dui 

dt 

due 

dt 

dnK+ 

dt 

dfiYfi 

dt 

dUKjCond 

dt 


+C nU2 - PneUi - 7^i,ads 

TC ^H2 "^Bjads "^gaSjTec T "^gaSjion T "^therm 

"^gaSjTec T "^gaSjion T "^Kjevap f + "^K+jColl 

T7^gas,rec "^gaspon T "^Kjevap (1 f+) "^KOjColl 

+’^K+,coll + ’^K0,coll — "^KiCvap 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


It is straightforward to conhrm that the total number of potassium atoms (neutral, ionized, 
and bound) is conserved in snch a model, i.e., nK,cond + ^^k+ +J^ko = ''^K.tot, where we assnme 
^K,tot = 3.04 X 10“^nH2 (Lodders 2003). We seek a steady-state solution by setting all time 
derivatives to zero. Note that the hve densities are, through the cross section modihcations 
J, all fnnctions of the grain charge Z, the sixth variable. The sixth eqnation needed to close 
the system corresponds to charge nentrality. If grains’ mean charge is +Ze, 


rig^Z nK+ + nyig+ - Ue = 0. (27) 

In practice, we vary Z till charge nentrality is achieved and the other densities are in steady 
state. 

More specihcally, we hnd the steady-state solntion for gas of density nH 2 and 
temperatnre T as follows. At a particular value of Z, we calculate the cross section 
modihcations for ions, electrons, and nentral atoms, using the approximation formulas 
of Draine & Sutin (1987). We assume the gas and dust temperatures are identical; in 
protoplanetary disks this is an excellent assnmption except for the nppermost < O.lgcm”^ 
snrface layers (Glassgold et al. 2004). With the coefficients known, we then solve for the 
species’ steady-state abnndances. Given guesses for the densities of K° and K+, we hnd the 
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electron density by first solving for Ui. 


Hi = 


C 


C 


f3 Tig TTgr r'^i^ads /3 Tig -f" TTgr TTOgj. C\ Si Ji 

This value is then substituted into the equation for the electron density, 

C^Ha T "^therm T ^2 ^Ha [/^ T k—2 ^Ha ^K+ T ^e,ads] ) 

to derive a quadratic equation for Ue- 

1/3 ("e.adi + k-1 WHa nK+)] Ke 


(28) 


(29) 


Up 


+ [^'i.ads («^e,ads + ^-2 ^^Ha ^K+) ” P (^^gr dTrOg^, j(T) + /Cs URa ^ko)] 

- t'i,ads [C ^^Ha + dvrOgj, j(T) + /Cs URa ^ro] = 0. (30) 

Once rip is found, rii is found using the above equation, as are the other densities, as follows: 

nK,tot {k2 URj + ^2x0,0011 /^) 


nK+ = 


k2n}i2 + k-2nn2np _|_ z/k+,co11 (1 — /'*') + ^ko.coIi/'*' + 

\ t^evap J V ^evap J ’ ' ' ’ i^evap 

(31) 


yields the K"*" density, in terms of which 

nK,tot — n-K+ ^1 + 


*^K+,coll 


UrO — 


1 + 


'^kO.coII 


and 


^K,cond 


nK,tot - nK+ 

\ '^k 0 ,co 11 / 


1 + 


(32) 


(33) 


3 'k 0 ,co 11 


In practice, at each guess for Z we iterate between solving for rip and solving for ni, uro 
and nR+, then updating rip. Performing 10 such iterations allows us to compute the net 
charge rii + nR+ + n^^-Z — rip for the given Z. At very large Z > 0 the net charge is positive, 
while at very large Z < 0 it is negative, so it goes to zero at some hnite Z. We use a 
bisection method with 80 steps to hud the value of Z that yields both charge neutrality and 
ionization equilibrium. 
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3. Results 

3.1. Canonical Case 

We have solved for the abundances of electrons, ions (Mg’*" and K"*"), as well as neutral 
and bound K atoms, using the chemical network described above. We start with our 
canonical case, characterized by nHj = lO^^cm”^, grains with a dust-to-gas mass ratio 0.01, 
uniform radius Ogr = 1 /rm, and work function W = 5.0 eV. Figure 2 shows the fractional 
abundances of electrons, Xg = Ue/nHa, as well as Xi, xk+, x^o, and a;K,cond, as functions of 
temperature. In cold gas, < 500 K, the ion and electron densities come from the balance 
between ionizations by energetic particles and collisions with dust grains. Because of their 
slower thermal velocities, ions take longer to reach the grains and are more abnndant, 
with Xi ~ 10“^®, compared to Xg ~ 10“^^ for electrons. The very low abnndances are a 
consequence of the low ionization rates associated with the radioactivities. Essentially 
all potassium is bound into solids, with few ions or nentral atoms in the gas phase. As 
the temperatnre increases, however, and especially K+ enter the gas phase. Some 
come directly from evaporation of K+ ions from grain snrfaces; because IP < W, most K 
atoms leave as ions. In addition, once the temperatures exceed about 750 K, electrons are 
produced effectively by thermionic emission. Above 1000 K, the potassium atoms evaporate 
from solids; most takes the form of neutral K° atoms in the gas-phase. Only in gas hotter 
than abont 1500 K do the K"*" ions compare in abundance to the nentral K° atoms. 



21 



600 800 1000 1200 1400 

T (K) 


Fig. 2.— Abundances of electrons (solid red curve), atomic (Mg) ions (dotted blue curve), 
K+ ions (dashed green curve), gas-phase K° atoms (dash-dot curve), and bound K atoms 
(dash-dot-dot-dot curve), as functions of temperature T, for the canonical case (gas density 
riHa = lO^^cm”^, dust-to-gas mass ratio 0.01, grain work function W = 5.0 eV. At low 
temperatures, non-thermal ionizations (from radioisotopes) dominate, but emission of K"*" 
ions from dust grains is signihcant for T > 550 K. Evaporation of K° atoms from grains, 
and their subsequent thermal ionization, contribute signihcantly to K"*" and electron densities 
for T > 750 K. For T > 1000 K, thermal ionizations and recombinations dominate and are 
nearly in thermodynamic equilibrium. For comparison, we also plot the electron abundance 
derived using the Saha equation (dashed red curve). The Saha equation overpredicts the 
ionization at temperatures ^ 800 K, and underpredicts the ionization at higher temperatures. 
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The bolded curve in Figure 3 shows the charge on dust grains as a function of 
temperature. Below about 700 K, when the ionization rate is low, grains typically hold 
only a few electron charges. The value of Z is such that charging by electron and ion 
collisions balance each other, including the effects of charge on the collision rates (i.e., Je, 
Ji). As the temperature increases, however, ion emission comes to dominate the charging, 
and the grains quickly become negatively charged; at 900 K the charge is about —700e. 
Because adsorption of gas-phase ions is unimportant, the grains charge until Wes is reduced 
enough that ion emission is suppressed and thermionic emission enhanced. In other words, 
ions are emitted until |Ze^/agrIF| approaches unity, which requires Z ~ 3500. At 900 K, 
\Ze^/a^r\ ~ 0.2 W. As the temperature increases, collisions of ions with grain surfaces begin 
to matter, and the grain charge is reduced, approaching the thermodynamic limit where 
|Ze^/agr| ~ kT, or Z ~ —100 (at 1500 K). 
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Fig. 3.— Charge on dust grains (in units of elementary charge e) as a function of temper¬ 
ature, for the canonical case but with various grain work functions. From top to bottom, 
W = 2.0, 2.5, ... 6.0eV. The bolded curve corresponds to W = 5.0eV. For low work 
functions (W ^ 4eV, thermionic emission dominates the charging, and dust grains are pos¬ 
itively charged; as temperature increases, thermionic emission becomes more effective, and 
grains become more positively charged. For large work functions {W > 4eV), ion emission 
dominates, and grains become negatively charged; as the temperature increases, collisions of 
K+ ions onto the grains becomes more effective, and grains become less negatively charged. 


The relative abundances of neutral and ionized K atoms also approach the 




24 


thermodynamic limit only at high temperatures. The bolded curve in Figure 4 shows what 
effective ionization potential energy E must exist so that the Saha equation would yield the 
correct ratio for nx+Ue/nKO, i.e., 

F' = — [In (nK+ Ue/nKo) — 35.42] kT. (34) 

At high temperatures it tends to approach the actual value of IP for potassium, 4.34 eV, but 
at lower temperatures E is lower, reflecting the fact that many K atoms evaporate as K"*" 
ions from grain surfaces, according to the Saha-Langmuir equation. At all temperatures, 
but especially below 800 K, emissions from dust grains generate more K"*" ions than the 
Saha equation would indicate. 
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Fig. 4.— Effective ionization potential, E, needed for the Saha equation to return the correct 
value of nK+ ne/nxo, as a function of temperature, for various values of the work function, 
W. The value of W (in eV) is used to label several of the curves. The case W = 5.0 eV is 
bolded. The horizontal dashed line refers to E = 4.34 eV, the value tht would obtain in the 
absence of species other than potassium. Because of emission of charges from dust grains 
(thermionic emission of electrons at low W, ion emission of K"*" at large W), the gas charge 
state is always greater than expected based on a pure Saha equation for potassium; that is, 
E < 4.34 eV. The effect is most pronounced and sensitive to W at low temperatures < 800 K. 
For T > 800 K, 3.5eV < E < 4.34 eV for W between 2 and 6 eV. For work functions close 
to the ionization potential of potassium (W ~ 4.34 eV), E 4.34 eV. 
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Figure 5 shows the timescales for the steady-state electron density to be affected 
by various ionization / recombination processes, as functions of temperature. For 
example, since IZthem is the rate per volume of nebula at which electrons are produced, 
hherm = /TZiherra- Likewise, the timescale for recombination of electrons with atomic 

ions is dehned as ne/{/3 nen{) = l/{/3ni). At low temperatures, T < TOOK, the dominant 
loss mechanism is adsorption of electrons onto dust grains, while the dominant production 
mechanism is ionizations of gas-phase molecules by energetic particles. These balance 
each other and the electron density turns over on timescales of roughly 10 seconds. At 
temperatures above 800 K, thermionic emission of electrons from grains dominates, and 
the electron density increases. Ionizations by energetic particles become less important 
compared to thermionic emission. But because the electron density is greater, the timescale 
for overturning the electron population rises, to above 1 hour at about 780 K, to about 
1 day at 870 K, and peaking at 10® s at 900 K. At higher temperatures still, the timescales 
for both adsorption and thermionic emission decrease, eventually falling below about 1 hour 
again for T > 1300 K. 
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Fig. 5.— Characterisic timescales for various processes affecting the electron abundances, as 
functions of temperature, for our canonical case. The dominant removal processes is always 
electrons colliding with grains (‘adsorption’, solid curve). At low temperatures (T < 800 K), 
ionization of H 2 by energetic particles (’ionize, H 2 ’, dashed curve) is the dominant charging 
process. At higher temperatures, emission of electrons from grains (’thermionic’, long dashed 
curve) dominates. At high temperatures the timescales for thermal ionization of K (dashed- 
dotted curve) and recombinations of K"*" ions and electrons (dotted curve) become shorter 
but never dominate. Recombination of atomic ions and electrons is slower than all other 
processes and is not plotted. Notably, for temperatures 800 K ^ T < 1300 K, the chemical 
timescales are longer than hours. 







3.2. Variation with Work Function 


We now investigate the effect of varying the work function of the solids making up 
the grains, keeping other parameters the same as our canonical case (nna = 10^"^cm“^, 
dust-to-gas ratio 10“^). Figure 3 shows the charge on each dust grain (in units of e) as a 
function of temperature T for various W. At low work function, W < 4eV, thermionic 
emission of electrons is very effective, whereas any potassium atoms leaving are likely to 
become neutral K° atoms. The grains are positively charged. At higher temperatures, 
the rate of thermionic emission is greater, and the dust grains are increasingly positively 
charged. For high work function, W > 4eV, thermionic emission is ineffective, but emission 
of K+ is effective. At higher temperatures, grains emit ions more effectively and so are more 
negatively charged. At temperatures ~ 800 to 1000 K, however, ion emission is as effective 
as it can be. At still higher temperatures, gas-phase ions’ collisions with dust grains are 
increasingly important, and the grains become less negatively charged. 

Figure 4 shows the effective ionization potential, E, needed for the Saha equation to 
yield our calculated value of nK+n,e/nK 0 , as a function of temperature, for various work 
functions W. At low temperatures this ratio depends sensitively on hF, as the ion and 
electron abundances are controlled by grain processes. Especially for high work functions, 
VF > 4 eV, ion emission dominates and increases nK+ above the levels predicted using the 
Saha equation. At low work functions, thermionic emission increases Ue- In gas hotter than 
about 800 K, thermal ionizations dominate and E is driven toward the ionization potential 
IP = 4.34 eV; however, this is only a trend, and the effective ionization potential falls in 
the range 3.5 eV < F < 4.34 eV. The departure of E from IP is most pronounced for W 
that are far from IP; for W ~ 4.34 eV, E 4.34 eV. Notably, though, E never appreciably 

exceeds IP, meaning the dust makes the gas better-ionized than would be predicted using 
the Saha equation alone. 
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Figure 6 shows the abundances of various species if the work function is increased to 
W = 6.0 eV. Compared to the case with W = 5.0 eV, the abundances as a function of 
temperature T are quite similar, with the notable differences that hko is much reduced, 
and nK+ rises with T much more steeply. This underscores that what is driving up nK+ 
for T < 800 K is ions evaporating from grains, and whether the K atoms leave as ions or 
neutral atoms is a strong function of W, or more precisely hF — IP. Because of the grains’ 
large negative charge, atomic ions’ abundances are suppressed. 
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Fig. 6.— Abundances of electrons (solid red curve), atomic (Mg) ions (dotted blue curve), K+ 
ions (dashed green curve), gas-phase K° atoms (dash-dot curve), and bound K atoms (dash- 
dot-dot-dot curve), as functions of temperature T, for the case with canonical parameters 
(see Figure 2), but with grain work function W = 6.0eV. The electron fraction predicted 
using the Saha equation is shown for comparison (dashed red curve). 


Figure 7 shows the abundances if the work function is decreased to W = 2.0 eV. Again, 
the general behavior with temperature is similar, except for T < 800 K. For the W = 2.0 eV 
case, production of K+ ions is considerably lower; according to the Saha-Langmuir equation, 
nearly all K atoms leave as neutral atoms. On the other hand, production of free electrons 
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by thermionic emission is much more effective. The electron fraction rises to ~ 10“^^ 
and plateaus there for temperatures between about 550 K and 900 K. Because of the large 
positive charge on each grain, the gas-phase ions’ abundances are elevated. 
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Fig. 7.— Abundances of electrons (solid red curve), atomic (Mg) ions (dotted blue curve), K+ 
ions (dashed green curve), gas-phase K° atoms (dash-dot curve), and bound K atoms (dash- 
dot-dot-dot curve), as functions of temperature T, for the case with canonical parameters 
(see Figure 2), but with grain work function W = 2.0eV. The electron fraction predicted 
using the Saha equation is shown for comparison (dotted red curve). 


Finally, Figure 8 shows just the electron abundance as a function of temperature, for 
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various values of the grain work function. At high temperatures, > 1000 K, the behavior of 
Xe with T is very similar, but the exact value does depend on the work function, although 
not monotonically with W. At lower temperatures the electron fraction is especially 
sensitive to the work function, being up to 8 orders of magnitude higher for low values 
{W ~ 2eV) than for high values {W ~ 6eV). 
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Fig. 8.— Electron fraction as a function of temperature, for various grain work functions. 
Other parameters are held at their values in the canonical case. The fiducial case {W = 
5.0 eV) is in bold. Electron fractions at low temperature (< 800 K) decrease as the work 
function is increased from 2.0eV to 3.0eV to 4.0eV (as labeled). As the work function is 
increased more, from 4 to 6 eV, the electron fraction tends to increase for temperatures 
T ^ 800 -900 K. 



3.3. Variation with Dust-to-Gas Ratio 


We now investigate the effect of varying the dust-to-gas ratio, keeping other parameters 
the same as our canonical case {W = 5eV, = lO^'^cm”^). Figure 9 shows the 
abundances of various species for the case of a dust-to-gas ratio 0.0001, a factor 100 below 
the canonical value. Figure 10 shows the abundances of various species for a dust-to-gas 
ratio of unity, a factor 100 higher than the canonical value. Comparing Figures 9 and 
2, for the reduced dust-to-gas ratio both the electron and ion densities are higher at 
low temperatures (T < 600 K), because both are produced by energetic particles and 
lost from the gas phase by adsorption onto grains. Comparing Figures 10 and 2, for the 
increased dust-to-gas ratio these densities are reduced below the canonical case. At high 
temperatures (T > 1000 K), all three dust-to-gas ratios lead to similar electron densities, 
and TT-e ~ nK+, because in that regime thermal ionizations of potassium dominate. At 
intermediate temperatures, there is a strong tendency for the electron and ion (both atomic 
ions and K+) densities to decrease with increasing dust-to-gas ratio, which makes sense 
given that they are lost mostly to collisions with dust grains. 
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Fig. 9.— Abundances of electrons (solid curve), atomic (Mg) ions (dotted curve), K+ ions 
(dashed curve), gas-phase K° atoms (dash-dot curve), and bound K atoms (dash-dot-dot-dot 
curve), as functions of temperature T, for the case with canonical parameters (see Figure 
2), but with dust-to-gas mass ratio 10“^. The electron fraction predicted using the Saha 
equation (dotted red curve) is shown for comparison. 
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Fig. 10.— Abundances of electrons (solid curve), atomic (Mg) ions (dotted curve), K+ ions 
(dashed curve), gas-phase K° atoms (dash-dot curve), and bound K atoms (dash-dot-dot-dot 
curve), as functions of temperature T, for the case with canonical parameters (see Figure 2), 
but with dust-to-gas mass ratio 1. The electron fraction predicted using the Saha equation 
(dotted red curve) is shown for comparison. 


Figure 11 shows the grain charge as a function of temperature, for various dust-to-gas 
ratios between 10“^ and unity. For all these cases, the charge on each grain is small 
at low temperatures, but becomes sharply negative as the temperature increases above 
some threshold; for the canonical ratio 10“^ that threshold temperature is 700 K. As 
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the temperature is increased further, the grain attains a maximum negative charge but 
then becomes less charged, eventually (above 1000 K in the canonical case) following 
Z ^ —500 (T/1200 K)“^ in all cases. For other dust-to-gas ratios the behaviors are the 
same but the thresholds occur at temperatures 100 K lower if the dust-to-gas ratio is 10“"^, 
or 200 K higher if the dust-to-gas ratio is unity. 

Finally, Figure 12 shows the electron fraction as a function of temperature for various 
values of the dust-to-gas ratio. As might be expected, is inversely proportional to this 
ratio at low temperatures (< 600 K), and independent of this ratio at high temperatures 
HOOK) where gas-phase reactions dominate both ionization and recombination. In 
between, the electron fraction is highly sensitive to the dust-to-gas ratio, being as much as 
12 orders of magnitude higher for small values (10“^), as compared to large values (unity). 
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Fig. 11.— Grain charge (in units of elementary charge e), as a function of temperature T, 
for various dust-to-gas ratios, from left to right 10“"^, 10“^, 10“^, 10“^, and 10°. The curve 
for the canonical dust-to-gas ratio, 10“^, is bolded. With increasing temperature, the grains 
hrst become more negatively charged as they emit K"*" ions, then become less negatively 
charged as they collide with ions. 
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Fig. 12.— Electron fraction vs. temperature for varying dust-to-gas ratios (as labeled). Other 
parameters are held at their values in the canonical case. At low temperatures (< 600 K), 
the electron fraction is inversely proportional to the dust-to-gas ratio. At high temperatures 
HOOK, the electrons come mostly from thermal ionization and their abundance is in¬ 
dependent of the dust-to-gas ratio. At intermediate temperatures, the electron fraction is 
highly sensitive (but still inversely related) to the dust-to-gas ratio. 
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3.4. Variation with Gas Density 

We now investigate the effects of varying the gas density, keeping other parameters 
the same as our canonical case {W = 5eV, dust-to-gas ratio 10“^). Figure 13 shows the 
electron fraction as a function of temperature for the case of a gas density nnj = lO^^cm”^, 
an order of magnitude below the canonical value. Figure 14 shows the electron fraction as 
a function of temperature for the case of a gas density nna = 10^® cm“^, a factor of 100 
above the canonical value. Figure 15 shows the average charge on grains as a function of 
temperature, for various gas densities. The same trends are apparent: at low temperatures 
(< 700 K), grains are nearly neutral, becoming very negatively charged at temperatures of 
about 800 to 1000 K, then becoming less negatively charged at higher temperatures. For 
higher gas densities, these transitions occur at slightly higher temperatures; e.g., a ten-fold 
increase in shifts the transitions higher in temperature by only about 60 K. Likewise, 
Figure 16 shows the electron fraction as a function of gas density. Again it is seen that the 
charging is relatively insensitive to the gas density. 
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Fig. 13.— Abundances of electrons (solid curve), atomic (Mg) ions (dotted curve), K+ ions 
(dashed curve), gas-phase K° atoms (dash-dot curve), and bound K atoms (dash-dot-dot-dot 
curve), as functions of temperature T, for the case with canonical parameters (see Figure 
2), but with gas density lO^^cm”^. The electron fraction predicted using the Saha equation 
(red dotted curve) is shown for comparison. 








42 



600 800 1000 1200 1400 

T (K) 


Fig. 14.— Abundances of electrons (solid curve), atomic (Mg) ions (dotted curve), K+ ions 
(dashed curve), gas-phase K° atoms (dash-dot curve), and bound K atoms (dash-dot-dot-dot 
curve), as functions of temperature T, for the case with canonical parameters (see Figure 2), 
but with gas density 10^® cm“^. 
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Fig. 15.— Charge on dust grains as a function of temperature, for our canonical values (see 
Figure 2), but for various values of rin^. 

4. Application to Disks 
4.1. Inner Edge of MRI Dead Zone 



Our chemical network allows us to calculate ionization levels in dusty plasmas at high 
temperatures. An important application is to determining the inner edge of the MRI dead 
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Fig. 16.— Electron fraction as a function of temperature, for our canonical values (see Figure 
2), but for various values of nHa- 

zone, that region of the disk where the MRI cannot lead to turbulence at the midplane. 
Jin (1996) and Gammie (1996) both realized that ionization levels at PPD midplanes are 
too low to permit good coupling of the gas to the magnetic helds, at least in regions where 
the temperatures are below about 1000 K. Closer to the Sun, at heliocentric distances 
inside roughly 0.1 to 1 AU, temperatures may exceed this threshold, and they argued that 
thermal ionization of alkali elements in the gas would lead to ionization levels high enough 
to permit the MRI. Because the mass flow changes abruptly at the dead zone edge, the 
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location of the transition bears significantly on the disk’s strnctnre and the interpretation 
of its spectral energy distribntion (e.g., Hasegawa & Pudritz 2010; Zhu et ah 2010). 


Using onr chemical network, we identify precisely what temperatnres are reqnired for 


snfficient conpling. In Fignres P^Tp^ we plot the Elsasser nnmbers associated with Ohmic 
dissipation (Aq), ambipolar diffnsion (Aa) and the Hall effect (Ah). The Elsasser nnmber 
v\l{'qfl) is the time in orbits for magnetic helds to diffuse across the fastest-growing MRI 
wavelength; va = H/( 47 rp)^A jg the Alfven speed, O is the Keplerian orbital frequency, and 
r] is the relevant diffusivity. The Ohmic, ambipolar and Hall diffusivities are straightforward 
to compute from the charged species’ abundances, gyrofrequencies around the magnetic 
held of strength B, and rates of collision with other species (Desch 2004; Wardle 2007). In 
Figures [T7]{^ we plot the three Elsasser numbers as functions of temperature T, for various 
combinations of heliocentric distance r, gas density ?7 ,h 2, and magnetic held strength, set by 
hxing the plasma beta parameter (3 = Pgas/(H^/87r). MRI turbulence is switched oh by an 
Ohmic Elsasser number Kq < 1 (Sano & Inutsuka 2001, Turner et ah 2007) or an ambipolar 
Elsasser number Aa < 1 (Bai & Stone 2011), while Hall Elsasser numbers Kh ^ 1 replace 
the turbulence with a laminar how in which magnetic torques nevertheless yield rapid 
angular momentum transport (Kunz & Lesur 2013; Lesur et al. 2014). 


We assume a 1 Mq star in what follows. Figure [T7| pertains to conditions at 0.1 AU. 
The top panel has nH 2 = 10^^ cm“^ and R ~ 1 G, appropriate for the disk atmosphere well 
above the midplane, and the bottom panel has nH 2 = 10^^ cm“^ and R ~ 1 G, appropriate 
for the midplane. As ionization levels increase with increasing temperature, the Elsasser 
numbers also generally increase. For the disk atmosphere, the last quantity to do so 
is clearly the one associated with ambipolar dihusion, which therefore is the dominant 
mechanism suppressing the MRI and setting the location of the dead zone inner edge. In 
the disk atmosphere at 0.1 AU, temperatures must exceed 1150 K for the MRI to operate. 
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At the midplane, ignoring the complications introduced by the Hall effect, ambipolar and 
Ohmic diffusion are equally important for setting the dead zone inner edge, which will lie 
where temperatures exceed about 950 K. 


Figure 18 pertains to conditions at 1 AU. The top panel has uh 2 = 10®cm“^ and 
B ~ 10“^ G, appropriate for the disk atmosphere, and the bottom panel has nH 2 = 10^^ cm“^ 
and B ~ 10“^ G, appropriate for the midplane. In the atmosphere, ambipolar diffusion is 
again seen to be the main process suppressing the MRI, and in the range of temperatures 
we consider, the MRI is never allowed. At the midplane, ambipolar diffusion once more 
dominates, and initiating the MRI requires temperatures above about 900 K. 


Finally, Figure 19 pertains to conditions at 0.1 AU, but during an epoch of increased 
density, such as might occur around the onset of an FU Orionis accretion outburst. The 
top panel has nH 2 = 10^^ cm“^ and R ~ 10 G, appropriate for the disk atmosphere, and 
the bottom panel has nH 2 = 10^® cm“^ and R ~ 10 G, appropriate for the midplane. In 
the atmosphere, ambipolar diffusion is again the dominant process suppressing the MRI; 
the Hall term is comparable. Temperatures must exceed about 950 K for the MRI to 
operate. At the midplane. Ohmic dissipation and the Hall effect are of similar strengths, 
and temperatures must exceed about 1070 K for the MRI to operate. 


If these examples are representative, then it would seem that the dead zone’s inner 
edge is set by the action of ambipolar diffusion, not Ohmic dissipation. In other respects, 
however, our results are similar to those that have been found using the Saha equation. In 
Figures [iTlfT^ we also plot the Ohmic Elsasser number under gas-phase collisional ionization 
of potassium, with the resistivity obtained from the standard formula r] = 

(Blaes & Balbus 1994; Equation 41). Our chemical network with thermionic and ion 
emission yields Ohmic resistivities lower, and Elsasser numbers higher, by an order of 
magnitude or more at temperatures near 1000 K. For the range of conditions we consider. 




initiating the MRI in disk atmospheres or midplanes generally requires temperatures of at 
least 900 to 1150 K, remarkably similar to the generally assumed value of about 1000 K. 


4.2. Stellar Magnetospheric Loading 

Young stars grow by accreting material from their disks. The gas flows to the star 
along the stellar magnetic held lines if the mass accretion rate is too low for the incoming 
gas to crush the stellar magnetosphere (Edwards et al. 1994; Bouvier et ah 2007). The gas 
could be loaded onto the stellar held either through the held lines’ dihusion into the disk 
(Miller & Stone 1997), through magnetic reconnection (Hirose et al. 1997), or through a 
Rayleigh-Taylor type instability (Romanova et al. 2008). Which scenario predominates, and 
how quickly the mass is loaded by each of the three mechanisms, depends on the magnetic 
dihusivity. We suggest that the transition from loading by dynamical instability to loading 
by magnetic dihusion occurs when the dihusivity exceeds where h is the disk density 
scale height. This is based on the Ending from ideal-MHD modeling that the plasma beta 
is about unity at the magnetosphere’s edge where loading occurs (Romanova et al. 2012). 
Consequently the scale on which magnetic and Coriolis forces balance, and fastest-growing 
vertical wavelength of the magnetorotational instability, will be about h. 

The mass loading takes place not too far from the co-rotation radius, where the 
Keplerian orbital period matches the star’s rotation period. Accreting T Tauri stars 
typically rotate about once a week (Bouvier et al. 1993), so the co-rotation radius around a 
1 Mq star is at r 0.07 AU. In Figure we present results near the co-rotation point for 
the Ohmic dihusivity rjo, Hall dihusivity rjn, and ambipolar dihusivity as functions of 
temperature. Three combinations of density and plasma beta parameter /3 = Pg/Pmag are 
shown: nm = lO^^cm”^, /3 = 1; nm = lO^^cm”^, /9 = 1; and rim = lO^^cm”^, /? = 100. 
We choose /3 = 1 for the hrst two panels, yielding a local held strength B that increases 
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slightly with temperature. In the third panel we select /3 = 100 since the stellar magnetic 
held at 0.07 AU is unlikely to exceed ten gauss, given a few kilogauss value at the stellar 
photosphere and a dipolar falloff with distance. The third panel’s parameters are thus more 
likely to be found in the higher accretion rate regime where the inhowing gas has enough 
ram pressure to force the stellar helds inward. 


Grey shading in Figure indicates diffusivities high enough that diffusive loading is 
likely to dominate. The transition from dynamical to diffusive loading takes place at a 
temperature in the range from 960 K (at nH 2 = lO^^cm”^) to 1080 K (at ?7 ,h 2 = lO^^cm”^). 
Ambipolar diffusion causes the transition except in the case with the highest density and 
weakest magnetic held, where Ohmic diffusion is stronger. The Hall term is also signhcant 
near the threshold in all three cases. Temperatures below the transition are quite plausible 
if accretion heating is weak. The thresholds are comparable to the 1060 K equilibrium 
temperature of a blackbody directly exposed to the starlight at 0.07 AU. Furthermore 
corotation lies beyond the dust sublimation front, so starlight can be absorbed before 
reaching 0.07 AU. In a passively heated protoplanetary disk, temperatures at 0.07 AU are 
predicted to be < 690K (Lesniak & Desch 2011). We conjecture that the stellar magnetic 
held may dihuse slowly into the hot gas nearer the star, but held lines that escape this 
region—either through reconnection or disruption of the disk’s inner edge by a dynamical 
instability—can quickly dihuse further into the disk and load gas from the disk on the 
stellar magnetosphere. 


Two points seem worth emphasizing. First, the close match between the dihusive 
loading threshold and the temperatures expected in passively heated disks means that the 
loading mechanism can switch depending on conditions. Second, at temperatures near 
the threshold, thermionic and ion emission reduce the dihusivity by orders of magnitude 
compared with the Saha equation. The grains’ work function thus likely impacts the 



character of diffusive loading on young stars’ magnetospheres. These effects should be 
included in future studies. 


4.3. Current Sheet Instabilities 

The treatment of high-temperature ionization in protoplanetary disks also can be 
applied to the “short-circuit” instability of current sheets. In protoplanetary disks with 
magnetorotational turbulence, it is reasonable to assume (and simulations show) that a 
commonly occurring magnetic geometry has one region of more-or-less uniform magnetic 
held adjacent to another region with the held lines oriented in the opposite direction (Sano 
et ah 2007; Hubbard et ah 2012). The surface separating the two is marked by steep 
magnetic gradients and high current densities, and is termed the current sheet. Ohmic 
dissipation can lead to heating in this sheet. Hubbard et al. (2012) argued that if the gas’ 
magnetic dihusivity rj decreases with increasing temperature then the heating can grow 
exponentially: increased temperature leads to lower dihusivity; the current concentrates 
in the lower-dihusivity layer, accompanied by more Ohmic heating; so the temperature 
increases further. The runaway heating can be quite rapid, potentially taking hours or 
less, depending on the ionization properties of the gas. Hubbard et al. (2012) and McNally 
et al. (2013) further argued that this sudden heating might be the mechanism behind 
the formation of chondrules, which are known to have been melted suddenly by some 
“hash-heating” process, then cooled and crystallized over the course of hours (Connolly 
et al. 2006). Whether current sheet instabilities can take place in protoplanetary disks, 
and whether they can grow fast enough to be the chondrule formation mechanism, are key 
questions that hinge on the ionization state of the gas at high temperatures. 

The growth rate of the short-circuit instability was calculated by Hubbard et al. (2012), 
who performed a linear analysis on a background state in which the magnetic held increases 
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linearly from —Bq/2 to +Bq/2 as x runs from —L/2 to +L/2, where x is the direction 
normal to the current sheet. They assumed that perturbations in the magnetic held, current 
density, etc., vary as exp (+Af + ikx), where k is the wavenumber of the perturbation and 
A is the growth rate. The growth rate is then found to satisfy the equation 

(Ool + ^diff + + ^diff (^cool ~ 

where fcooi is the timescale over which emission of radiation cools the region; we specify 
its value below. Here tdiff = l/(^^ho) is the timescale for the magnetic held to dihuse 
one wavelength, and ? 7 o is the magnetic dihusivity in the background state. Other zero 
subscripts likewise indicate conditions in the background state. The timescale can be 
dehned as 


t 


V ~ 


d\nri\ ^ 


(36) 


where 
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(37) 
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is the timescale for Ohmic dissipation to heat the gas signihcantly. We have written fheat in 

terms of the plasma beta parameter in the zero-order state, (Aq = Po/{BQ/8n), which makes 
clear that the heating rate is directly tied to the rate at which magnetic held lines dihuse 
into the current sheet. In the usual case where r] decreases with increasing temperature, 

B > 0. In that case, unstable growth is assured if B < tcooi- 


To determine fcooi we imagine a slab of thickness L, heated to temperature T by 
Ohmic dissipation. The slab will cool as it emits radiation that must dihuse through the 
surrounding gas. If this gas has density p, heat capacity Op, and opacity (per mass of gas) 
K, the appropriate timescale is the radiative dihusion time, 

3 p^ Op K 2 


frd 


(38) 


64 TT^ aT^ ’ 

(Mihalas & Mihalas 1984; Kippenhahn & Weigert 1990; Desch & Connolly 2002). So that 
the scenario is relevant to chondrule formation, we adopt p = 1 x 10“®gcm“^, T = 1500 K, 
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and C-p = 8 X lO^ergg”^ K“^. As for the opacity, a reasonable valne wonld be 26cni^g“^ 
if dust were present (at T = 1500 K), but since micron-sized dust generally vaporizes 
within minutes in the chondrule-forming region, a more reasonable opacity would be that of 
chondrules alone, k ^ 0.03 cm^ g“^ (Morris & Desch 2010). We hnd Ad = 5200 (L/10^ km)^ s 
if K = 26cm^g“^, and Ad = 6(L/10^km)^s li k = 0.03cm^g“h Note that McNally et ah 
(2013) took densities p ~ 10“^gcm“^, which would make cooling slower by 4 orders of 
magnitude; but even in a disk with 10 times more material than the minimum mass solar 
nebula (Weidenschilling 1977), these densities are only achieved far inside 1 AU, so some 
outward transport process would be needed to place the chondrules in the asteroid belt. 

At the lower densities that seem more likely to prevail, the current sheet cools by radiative 
diffusion in a fcooi of seconds to hours. 


With the cooling timescale in hand, it is simple to hnd the ratio of and tcooh 


tr. 


f dlnr]\ 


-1 
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647r^ aT^ 


(39) 
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Note that the ratio does not depend on the size of the slab: the heating timescale relies on 
diffusion of magnetic hux, and the cooling timescale relies on diffusion of radiation, so both 
scale as L^. Unstable modes require this ratio to be < 1, or 


V > hmin = 


dlnr] 


d\nT 
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Po 647r^ aT^ 


(40) 


2(7 — 1) 3 p^ Cp K 

For reasonable parameters {p = 10“®gcm“^, T = 1500 K, 7 = 7/5, Ap = 8x lO'^ergg”^ K“^), 
and assuming Pq = 10 and k, = 26cm^g“^, r/r nin ~ 2.4 x 10^'^ |(91n?7/(91nT|~^ cm^s“^. (or 
~ 10^ times higher if dust evaporates). The resistivity must be high enough that the current 
sheet can be continually fed with magnetic held lines by dihusion, faster than radiation can 
cool the gas. We can also express this requirement as a maximum electron fraction Xe, using 
the following relationship: 

rp \ 1/2 


7] = 9060 


1500 K 


^ cm^ s ^, 


(41) 
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(Hayashi 1981; Balbus & Terquem 2001). We find the unstable modes require 
Xe < a;e,max ~ 3.7 X 10“^^ |(9lu?7/(9luT|, or three orders of magnitude smaller if dust 
evaporates. 

For our canonical case depicted in Figure 2, if T < 750 K there is no variation in the 
electron density with temperature, in which case Xe^max ~ 0, and no unstable modes can 
be expected. For T > 750 K, |(91n?7/(91nT| ^ 100 (the electron density rises 8 orders of 
magnitude as T increases from 750 K to 900 K), so Xe^max ~ 3.7 x 10“^^ (or three orders 
of magnitude smaller if dust evaporates). This means the instability can proceed until 
T ^ 900 K. For hotter temperatures the electron density is too high for held lines to diffuse 
into the current sheet fast enough to overcome the radiative cooling. So, for the short-circuit 
instability to act, the ambient temperature of gas surrounding the current sheet must fall in 
the narrow range 750 K < T < 900 K. The case with grains with work function W = & eV is 
similar to the canonical case. For the W = 2eV case, |(91n?7/(91nT| ^ 14 for T between 600 
and 800 K. We then hud Xe^max ~ 5 x 10“^^ or so. A wider range of temperatures, extending 
down to 650 K, is then allowed to be unstable. Again, if dust evaporates in the chondrule 
formation region, then the faster radiative cooling means the maximum allowable electron 
density is much lower, by about three orders of magnitude. 

A second requirement for the current sheet instability is that a current sheet can be set 
up in the hrst place. This in turn requires the existence of magnetic turbulence in the disk, 
almost certainly from the magnetorotational instability (MRI). The MRI is suppressed if 
the Ohmic Elsasser number Aq = {v\/Q)/r] is less than unity, or if > v\/Q, where is 
the Alfven number and O the Keplerian orbital frequency (Jin 1996). From the definition of 
the plasma beta number, v\ = (2//d) where C is the sound speed; in disks in which the 
MRI has saturated, typically (3 ~ 10^ (Miller & Stone 2000; Hirose et ah 2006). Therefore 
the existence of the MRI requires rj < 0.02 0^/0. Note that the existence of a current sheet 
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demands a low value of the magnetic diffusivity, while the short-circuit instability requires 
a high value of 77 . For reasonable parameters (say, T ^ 800 K, at 1 AU in the disk), the 
MRI requires rj < 2.8 x lO^^cm^s”^, or Xg ^ 2.5 x 10“^^. For our canonical case, there 
is an exceedingly narrow range of electron fractions in the ambient nebula that can both 
lead to current sheets and the short-circuit instability: 2.5 x 10“^^ ^ Xg ^ 3.8 x 10“^^, 
corresponding to a very narrow temperature range about 10 K wide, centered on about 
900 K. As the work function of solids is reduced, the range narrows. For the W = 2eV 
case, there is in fact no overlap; the electron fractions needed for runaway growth are 
not compatible with the operation of the MRI and the existence of current sheets. If 
dust evaporates in the chondrule formation region, there is no overlap in conditions for 
any value of W. The net effect of all these factors is that the gas ionization needed to 
generate magnetic turbulence and to initiate the short-circuit instability are unlikely to be 
simultaneously met in protoplanetary disks. 

Finally, we point out one additional problem for the short-circuit instability, which 
is its need to grow faster than the chemical timescales will allow. The timescales in the 
problem must be < tcooi, which we argue is no more than 5000 s (for p = 10“®gcm“^ 
and no evaporation of dust). The instability requires the gas resistivity to decrease as the 
temperature increases, on these timescales. The steady-state electron abundance indeed 
increases with temperature, but the ionization chemistry takes time to reach steady state. 
According to Figures 2 and 5, the electron abundance starts to rise with temperature above 
750 K or so, but the chemical timescales to achieve these higher electron abundances also 
rise. Above about 800 K, the chemical timescales exceed 5000 s and approach 10^ s. This 
effectively reduces |(91n77/(91nT|, potentially by 1 to 2 orders of magnitude, making the 
instability more difficult to initiate and grow. 


In summary, the short-circuit instability requires a rise in electron abundance with 
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temperature, together with an initially large magnetic diffusivity so that magnetic heating is 
faster than radiative cooling. However, generating cnrrent sheets in the first place reqnires 
a small magnetic diffusivity. Depending on snch parameters as gas density, opacity, grain 
work fnnction, etc., there is either no way to simnltaneously satisfy these reqnirements, 
or they are satished only in a narrow (perhaps 10 K) temperatnre range, typically above 
700 K. The short-circnit instability therefore appears to reqnire rather special conditions. 
Moreover, the reqnirement that the temperatnre start above 700 K is inconsistent with the 
presence of primary sulfnr in chondrnles, which indicates ambient temperatnres <650 K 
(Lodders 2003). Even if short-circnit instabilities can arise from current sheets, they do not 
seem compatible with the chondrnle formation mechanism. 


5. Conclusions 

The ionization state of the gas hotter than about 500 K determines how well 
protostellar disks conple to their internal magnetic helds and to the helds of the central 
star. Properly calculating the ionization state is especially important for assessing where 
the magnetorotational instability can operate, which bears on the qnestion of how the 
disks evolve over time. It bears on whether disk gas can be loaded onto a star’s magnetic 
held. It also bears on whether current sheets can dissipate magnetic energy via the cnrrent 
sheet instability. Previous treatments of ionization at high temperatnres have assnmed that 
potassinm atoms are in the gas phase and that they are ionized according to the Saha 
eqnation. This approach ignores the role of dust and the fact that ions and electrons strike 
grains far more often than they meet one another in the gas and recombine. Fnrthermore 
gas-phase collisions are not the only ionization process. At high temperatnres, the grains 
also emit electrons and ions into the gas. Balancing thermionic (and ion) emission against 
adsorption of electrons (and ions) on the grains yields the Saha-Langmnir eqnation, which 
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is similar in form to the Saha equation, but instead of the first ionization potential of 
the gas-phase atoms (e.g. 4.34 eV for potassium) involves the work function of the grain 
material. Work functions typical of the solid materials in protoplanetary disks are near 
5 eV (silicates at 5.0 to 5.4 eV, FeNi metal at 4.4 eV, graphite 4.62 eV). These exceed 
the hrst ionization potential of potassium, but because the grains present so much surface 
area, the disk gas is more ionized (with higher electron density) than expected based on 
the Saha equation. These processes are important at temperatures between 700 K and the 
dust sublimation threshold. We are the hrst to identify their role in ionizing the gas in the 
central regions of protoplanetary disks. 

We have incorporated thermionic and ion emission into a chemical network and solved 
for the steady-state abundances of electrons and other species as functions of temperature 
and density. Unlike previous analyses, we do not automatically assume that potassium 
atoms reside in the gas phase; in fact. Fodders (2003) predicts a 50% condensation 
temperature of 1006 K, meaning that below about 1000 K, most potassium is locked up 
in solids. We mimic this effect by using an evaporation energy 3eV that potassium 
atoms must acquire to leave the grain surface. Thus above 800 K, at least one in 10"^ 
or 10® potassium atoms is in the gas phase rather than locked up in the grains. This 
treatment is approximate, but even if all the sodium were hrmly bound to the grains at 
lower temperatures, cesium atoms are in the gas phase above 850 K (Fodders 2003) and 
their abundance is only 10*^ times less than potassium. Some alkali atoms thus occur in the 
gas below 1000 K. 

The most signhcant result from our canonical case is a major rise in the electron 
fraction at temperatures between 700 K and 900 K. We hud that qualitatively, the 
run of electron fraction with temperature is relatively insensitive to the work function 
(Figure 8), gas density (Figure 16), and dust-to-gas mass ratio (Figure 12). Quantitatively, 
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the temperature at which the electron fraction rises sharply does vary somewhat: if the 
dust-to-gas mass ratio is increased (decreased) two orders of magnitude, the rise in electron 
density is shifted about 100 K to higher (lower) temperatures (Figure 12), with similar 
results expected for changes in grain size. Likewise, increasing (decreasing) the gas density 
by two orders of mangitude causes the rise in electron density to occur at temperatures 
about 100 K higher (lower) (Figure 16). A notable feature of the results is that the 
electron density can almost be computed using the Saha eqnation for potassinm, but one 
must use an ionization potential less than 4.34 eV. Above 800 K, for plausible grain work 
functions, values between 3.5 and 4.0 eV yield the correct ionization state, and that state 
is signihcantly more ionized than at lower temperatnres. As a robust result, then, thermal 
ionization can be said to be signihcant above about 800 K. 

A notable feature of the high-temperature regime is that the grains become rather 
highly charged through directly emitting electrons (or ions). The grain charge can be 
positive or negative, bnt for the most likely work function values near 5 eV, the grains emit 
ions and are left with negative charges. The charging is strong above 800 K, and grains of 
radius 1 pm reach charges of about —700 e at 900 K (Figure 3). While much larger than the 
grain charges found neglecting thermionic and ion emission, this is eqnivalent to an electric 
potential of just 1 volt. The charging is strong enough to affect the grains’ coagulation at 
high temperatnres. The Coulomb energies that must be overcome for grains to collide are 
~ 10^ eV, or ~ 10“®erg, so micron-sized grains coming together slower than 10cms“^ will 
be mntnally repelled by the electrical force. The effects of grain charging on the coagulation 
of grains have been investigated by Okuzumi (2009). We also point ont that if there were 
any grains composed of materials with work fnnctions < 3.5 eV, then their positive charges 
would enhance their collision rates with the negatively charged grains. 


We have applied our charging model to the qnestion of the MRI dead zone’s inner 
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edge. At low temperatures the gas is insufficiently ionized to permit the coupling of the 
gas to the magnetic held that is fundamental to the MRI. We hnd that for reasonable 
disk parameters (gas density, magnetic held strength), the temperature needed for the 
MRI to operate is set by ambipolar rather than Ohmic dihusion. The typical threshold 
temperatures are in the range 900 to 1100 K. At the high densities found near the midplane, 
the thresholds are comparable to or up to 100 K less than when neglecting the thermionic 
and ion emission and considering only the Ohmic resistivity. At the lower, atmospheric 
densities, ambipolar dihusion makes the thresholds up to 300 K higher than those from the 
Saha equation and Ohmic resistivity. Nevertheless, the results validate the common belief 
that thermal ionization requires temperatures above 1000 K. We likewise have investigated 
the temperatures required to dihusively load the disk mass onto the stellar magnetosphere. 
We hnd once again that temperatures below 1000 K are required, and that the threshold is 
set by ambipolar dihusion. 

The high-temperature ionization also has implications for the onset of the “short-circuit 
instability” and its role in chondrule formation. This instability relies on the resistivity of the 
gas decreasing with increasing temperature, so that dissipation of magnetic energy during 
reconnection events can accelerate. We conhrm that above a critical threshold temperature 
(about 800 K in our canonical case, but potentially as low as 650 K), this condition is met. 
On the other hand, the presence of primary sulfur (Rubin et ah 1999) in chondrules suggests 
that they condensed in regions cooler than 650 K, the condensation temperatnre of sulfnr. 
In addition, we show that there are two extra and seemingly incompatible constraints with 
which the short-circnit instability must comply. The instability requires the resistivity to 
change on timescales faster than the gas can cool. This in turn imposes a lower limit on 
the resistivity, and an upper limit on the electron fraction, so that field lines can diffuse 
rapidly enough into the reconnection region; typically Xe must not exceed abont 10“^^. The 
short-circuit instability also reqnires magnetic tnrbulence to set np the current sheets of 
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oppositely oriented field lines. Assuming magnetic turbulence is generated by the MRI, this 
in turn imposes an upper limit on the resistivity, and a lower limit to the electron fraction; 
typically must be 10“^^ or greater. For typical conditions, only a small range of electron 
fractions satisfy both conditions simultaneously, while for many parameters (especially if 
dust evaporates in the chondrule-forming region), no conditions simultaneously lead to the 
formation of current sheets and allow the short-circuit instability. This and the constraint 
that chondrule precursors started at temperatures < 650 K to retain primary sulfur would 
seem to imply that current sheets and the short-circuit instability are probably not the 
chondrule formation mechanism. 

We note that thermionic and ion emission may be important to the ionization state 
of gas in the atmospheres of exoplanets, especially “hot Jupiters”. Ohmic dissipation of 
magnetic energy in their atmospheres has been invoked as the explanation for the inflated 
radii (above model predictions) of hot Jupiters such as HD209458b and HD189733b 
(Batygin & Stevenson 2010; Perna et al. 2010) and may affect atmospheric circulation 
(Menou 2012; Rogers & Komacek 2014). At present, the resistivity of the gas in such 
models is computed essentially using equation 41 (as in Balbus & Terquem 2001), with 
the electron density found using the Saha equation for various gas-phase atoms (Batygin 
& Stevenson 2010; Rogers & Komacek 2014). However, hot Jupiters, HD209458b and 
HD189733b in particular, are predicted to have cloud layers containing condensates such as 
TiO (Fortney et al. 2010). TiO has a work function of about 4.5 eV (Fomenko 1966; Greiner 
et al. 2012) and a cosmochemical abundance equivalent to a dust-to-gas ratio ~ 10“® 
(Fodders 2003). At pressures of 1 bar, and temperatures of 1000 K, the gas density would 
be UHa ~ 6 X 10^®cm“^. We therefore expect that thermionic and ion emission from TiO 
particles, and adsorption of electrons and ions onto them, will control the ionization state 
of the gas, which will then not be described using a Saha equation, above some threshold 
temperature. Comparing to Figures 12 and 16, we expect that the high density will raise 
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the threshold temperature, while the low dust-to-gas ratio would lower the threshold 
temperature, but that thermionic emission would dominate the charging of the gas above 
about 800 to 1000 K. Naively, we expect thermionic emission will lead to higher electron 
fractions and lower Ohmic dissipation than use of the Saha equation alone, but exoplanet 
atmospheres are far enough out of the regime considered in this paper that other details 
may play a role. In any case, we recommend that this charging mechanism be considered 
when calculating Ohmic dissipation in the atmospheres of hot Jupiters. 

We also recommend that future effort be directed to measuring the work functions 
of astrophysically relevant materials, especially those found in protoplanetary disks and 
exoplanet atmospheres. We culled information from a variety of decades-old sources, but a 
dedicated effort should be made to consistently compare thermionic and ion emission from 
silicates, FeNi, graphite particles, and other astrophysical solids. 

Where solids exist, they can catalyze the recombinations of electrons and ions on 
their surfaces, rendering the Saha equation incorrect. But they also emit electrons or ions 
directly into the gas, at rates that depend on their material properties. In the end, the 
ionization fraction increases faster with temperature than under the Saha equation. For 
relevant conditions, we hnd that the ionization rises signihcantly above 800 K. 

We thank the organizers of the August 2014 meeting “Non-ideal magnetohydrody¬ 
namics, stability, and dissipation in protoplanetary disks” at the Niels Bohr International 
Academy in Copenhagen, for providing a forum for our collaboration. This research was 
carried out in part at the Jet Propulsion Laboratory, California Institute of Technology, 
under a contract with the National Aeronautics and Space Administration, and with the 
support of grant number 13-OSS13-0114 from the Solar Systems Origins program. 
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Fig. 17.— Ohmic Elsasser number Aq (red solid curve), ambipolar Elsasser number Aa (blue 
dashed curve), and Hall Elsasser number Ah (green dotted curve) as functions of temperature 
at 0.1 AU from a Solar-mass star. The red dotted curve shows the Ohmic Elsasser number 
computed from the Saha equation for potassium, neglecting thermionic and ion emission. 
Different values of the density and plasma beta parameter in the top and bottom panels 
(top: nH2 = 10^^ cm“^, /3 = 1; (bottom: nH2 = 10^"^cm“^, (3 = 100) correspond to the disk 
atmosphere and midplane, respectively. 
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Fig. 18.— Same as Figure 17, but for heliocentric distance 1 AU. Ohmic Elsasser number 
Aq (red solid curve), ambipolar Elsasser number Aa (blue dashed curve), and Hall Elsasser 
number Ah (green dotted curve) are plotted vs. temperature T. The red dotted curve 
shows the Ohmic Elsasser number computed from the Saha equation for potassium. The 
top and bottom panels correspond to the disk atmosphere and midplane, respectively (top: 
UH2 = 10®cm“^, /S = 1; bottom: UH2 = 10^^ cm“^, (3 = 100). 
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Fig. 19.— Same as Figure 17, except that the densities and plasma beta parameters cor¬ 
respond to the disk atmosphere and midplane at the onset of an accretion outburst, when 
gas and dust have piled up near 0.1 AU (top panel: nm. = lO^^cm”^, /3 = 1; bottom: 
?7 ,h2 = 10^®cm“^, = 100). Ohmic Elsasser number Aq (red solid curve), ambipolar El- 
sasser number Aa (blue dashed curve), and Hall Elsasser number Ah (green dotted curve) 
are plotted vs. temperature T. The red dotted curve shows the Ohmic Elsasser number 
computed from the Saha equation for potassium. 
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Fig. 20.— Magnetic diffusivities vs. temperature T, for fiducial grain parameters and three 
different combinations of gas density and plasma beta parameter. In each panel we plot the 
Ohmic (solid red curve), ambipolar (dashed blue curve), and Hall diffusivities (dotted green 
curve; the Hall diffusivity is negative on the thicker portion). Electron and ion abundances 
are derived using our chemical network including thermionic emission. Also plotted (red 
dotted curve) is the Ohmic diffusivity derived from the Saha equation for potassium. Gray 
shading indicates diffusivities greater than high enough that disk gas ought to be loaded 
onto the stellar helds primarily by diffusion. The threshold temperature is about 1000 K 
and is set by ambipolar diffusion except in the highest density case. 
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